Observation of the dynamic Jahn- Teller effect in the excited states of 
nitrogen-vacancy centers in diamond 
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The optical transition linewidth and emission polarization of single nitrogen-vacancy (NV) centers 
are measured from 5 K to room temperature. Inter-excited state population relaxation is shown 
to broaden the zero-phonon line and both the relaxation and linewidth are found to follow a 
dependence for T < 100 K. This dependence indicates that the dynamic Jahn- Teller effect is the 
dominant dephasing mechanism for the NV optical transitions at low temperatures. 
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The negatively-charged nitrogen-vacancy (NV) center 
in diamond has attracted much scientific interest due 
to its unique spin and optical properties. Long coher- 
ence lifetimes of the ground electron spin state at room 
temperature, up to milliseconds in recent reports [T], 
have enabled fundamental studies of coherent electron- 
electron f5] and electron-nuclear coupling [31 IH [S] of a 
small number of spins in a solid. Spin-selective optical 
transitions enable readout of the electron spin at room 
temperature which is critical for both high sensitivity 
magnetometry [5] and quantum information processing. 
Moreover, at cryogenic temperatures, it is proposed that 
these spin-dependent optical transitions [7] H] could pro- 
vide an interface between spins and photons as needed 
in schemes for scalable quantum computation [9l [10] and 
quantum communication |lll I12j . The success of such 
schemes depends on the coherence properties of the emit- 
ted photons which in turn are determined by the coher- 
ence properties of the NV excited states. 

Dephasing of the excited state of an optical transition 
manifests itself in an energy broadening of the transi- 
tion. For solid-state defects, the broadening can be ob- 
served in the zero-phonon line (ZPL), the optical transi- 
tion in which no net phonon is emitted or absorbed. Typ- 
ically the ZPL width exhibits a dependence on tem- 
perature T, and the dephasing mechanism is attributed 
to electron-phonon scattering mediated by a quadratic 
electron-phonon interaction |131 [H] . In this work we first 
measured the polarization of the emitted ZPL photons 
as a function of T. The NV excited state consists of 
an orbital doublet and a change in the emitted polar- 
ization indicates population transfer between the two or- 
bital states. We found that the emission polarization is 
strongly temperature dependent and that the population 
relaxation rate between the orbital states varies as . 
Next we measured the temperature dependence of the 
ZPL width. A dependence was again found, indicat- 
ing that for T < 100 K the ZPL broadening is mainly 
determined by excited state population transfer. Most 
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FIG. 1: (a) Electronic energy level diagram for the NV cen- 
ter in the presence of a strain field, (b ) PL intensity of the 
Ex and Ey peaks as a function of half waveplate (HWP) an- 
gle for NV3. (c) ZPL spectra of a single NV excited with 
Y-polarization. black: Collecting both polarizations, red: 
Collecting Y polarization, blue: Collecting X polarization. 



significantly, the dependence observed is evidence of 
the dynamic Jahn-Teller (DJT) effect in the NV excited 
states, as explained below. 

The NV center in diamond has trigonal symmetry 
(point group C3v) and consists of a single substitu- 
tional nitrogen with a nearest neighbor carbon vacancy. 
The NV axis {z axis) can point along any of the four 
(111) crystallographic axes. The ground state has 
•^^2 symmetry (spin-triplet/orbital-singlet) and is split 
into an Sx, Sy doublet 2.87 GHz above an Sz singlet. 
These are connected by optical transitions to the spin- 
triplet /orbital-doublet '^E excited state. The excited- 
state structure involves spin-orbit and spin-spin interac- 
tions as well as a linear strain splitting of the two orbital 
states [H]. The diagram in Fig. [T^ applies to the case 
where strain dominates over these other interactions such 
that the two orbital branches and Ey are well sepa- 
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rated in energy. Here, x and y are mutually orthogonal 
axes in a plane perpendicular to the NV axis with an 
angle determined by the strain tensor [16] . The ZPL op- 
tical transitions to and Ey have orthogonal, linear 
polarization selection rules [T7] and are labeled X and 
Y. Excitation along the NV axis with X-polarized light 
should therefore produce X-polarized photoluminescence 
(PL), unless either population relaxes between the states, 
or else the selection rules are degraded, as occurs for ex- 
citation or collection through the phonon sidebands. 

In the first experiment we measured the ZPL emission 
polarization as a function of temperature. Five centers 
were studied in a (lll)-oriented type Ila natural dia- 
mond sample chosen for its low NV density. The stud- 
ied centers had their NV {z) axes oriented parallel to 
the excitation and collection path, as determined from 
their polarization anisotropy [18j[19]. For NVl-2, no ex- 
ternal stress was applied, but random strain fields were 
present due to local defects in the crystal. For NV3- 
5, external stress was applied by mounting the sample 
between copper plates which contracted during the cryo- 
genic cooldown process. All experiments were performed 
in a coldfinger vacuum cryostat. 

The polarization visibility V = (/x — Iy) / (Ix + ly), in 
which /x(Y) is the intensity of the X(Y) polarized emis- 
sion from the '^E state, was measured for both X and Y 
polarized 532 nm excitation. Alignment of linearly po- 
larized excitation light to the principal axes of the NV 
center was performed by using a 'pump' half waveplate 
directly before the final imaging objective. To determine 
the proper waveplate angle, the X/Y transition intensity 
was measured as a function of waveplate angle as shown 
in Fig. [ijj. The intensity extrema occur in this curve 
when the excitation field is either X or Y-polarized. In 
a system with ideal polarization selection rules without 
excited-state relaxation, V at these extrema should be 
equal to one, but we measure |y| < 0.5 in all cases. A 
polarizer inserted before the spectrometer confirmed that 
the emission lines are in fact well polarized as shown 
in Fig. [TJ;. This, combined with the modest visibility 
temperature dependence at 10 K (see below), indicates 
that the main factor limiting the visibility in Fig. [l]3 is 
imperfect excitation selection rules when exciting non- 
rcsonantly through the high energy phonon sidebands. 

The polarization visibility, plotted in Fig. [2^, shows 
a strong temperature dependence from 10 — 40 K indi- 
cating that in this range the polarization decay rate is 
comparable to the radiative decay rate. A fit of the low 
strain NVl-2 polarization data to a simple three level 
model is shown in Fig. [2^. In this model, depicted in 
Fig. [2]3, the two excited states are spht by energy A and 
the radiative decay rate is r = (12.5 ns)~^ [20] [21]. The 
model includes the excitation rate n to the selected ex- 
cited state plus a rate a x n to the other excited state. 
The relaxation rate b from states 3^2 is taken to have a 
temperature dependence b/r = ciT^ . The temperature 




FIG. 2: (a) Temperature dependence of the polarization visi- 
bility V = (/x-/y)/(/x-I-/y) for 5 NVs. The solid black curve 
corresponds to a fit of NVl-2 data to the three-level model 
depicted in (b). For the spectrally resolved NV3-5, the deter- 
mination of the strain axis is discussed in the text. For NVl-2, 
an 'analyzer' half waveplate followed by a fixed polarizer was 
inserted before the spectrometer to measure both components 
of the spectrally unresolved X/Y ZPL peaks. The stain axis 
was determined by measuring the X/Y intensities obtained 
from rotating the 'analyzer' waveplate for specific angles of 
the 'pump' waveplate. Strain splittings for NVl-5 were 8, 9, 
44, 54 and 81 GHz respectively, (b) Three level model used 
to fit the temperature dependent polarization relaxation for 
Y-polarized excitation. For X-polarized excitation, a and an 
are reversed. 



dependence for each NV is fitted using parameters ci , N 
and a, and good agreement with experiment is obtained 
yielding ci = (5.1 ± 4.2) x 10"^ K"^, iV = 5.0 ± 0.3 and 
a = 0.40 ±0.01 [ig. 

In the second experiment we measured the linewidth 
of the ZPL through photoluminescence excitation (PLE) 
spectroscopy. In this measurement, photoluminescence 
into the phonon sidebands was detected as a tunable 
external-cavity diode laser operating at 637 nm was 
scanned across the ZPL resonances. Before each scan, 
a 532 nm excitation pulse was applied to reverse photo- 
ionization which eventually occurs with 637 nm excita- 
tion alone. These measurements were performed on a 
commercially-obtained (lOO)-oriented synthetic diamond 
sample (Element 6, Electronic Grade CVD) chosen for its 
high purity and relatively small spectral diffusion. Two 
NV centers were studied, labeled NV 6 and 7. 

At the lowest temperatures, with a single excitation 
frequency, typically only a single PLE line is seen corre- 
sponding to the vrLs = 0, Ex transition, while all other 
transitions arc hidden due to optical pumping [71 151 [^. 
Transitions to the lower orbital branch {Ey) can be re- 
vealed by applying 2.9 GHz modulation to the excitation 
laser, as shown in the bottom plot of Fig. [3^. The peaks 
at ~ 10 ± 2.9 GHz correspond to the rus = 0, E^ transi- 
tion excited by the laser and modulation sidebands, while 
the peaks at ~ 17 and 20 GHz represent two situations 
in which the laser and modulation sidebands can simul- 
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FIG. 4: Log-log plot showing the corrected PLE linewidths 
vs. temperature from NV6 and NV7, as well as PL linewidths 
from NV3-5 (PL). The solid line is the fit described in 
the text. The strain splitting for NV7 was 7.2 GHz and was 
undetermined for NV6. 



FIG. 3: (a) PLE spectra from NV7 at several temperatures. 
Transitions are labeled according to the excited orbital state, 
(b) Individual PLE scans of NV7 illustrating spectral diffu- 
sion. Intensity scale: (white) to 8600 s~^ (black), (c) Sum 
of scans (blue) using the procedure described in the text to 
remove spectral diffusion, and a Lorentzian fit (red). 

taneously drive vris = and m^. = ±1, Ey transitions. 

At the lowest temperatures, linewidth estimation is 
complicated by spectral jumps and blinking of the PLE 
resonance, illustrated in Fig. [3]d. To remove this ef- 
fect from the calculated linewidth, individual scans were 
shifted according to their first moment before summing. 
A combined spectrum obtained in this way is shown 
in Fig. [3];. The raw linewidth was then obtained from 
a Lorentzian fit, also shown. To avoid photobleaching 
and power broadening we used weak excitation powers: 
for example for NV 7 at low temperature the excitation 
power was 17 nW (focused to a 1 /im spot), while the sat- 
uration power was estimated to be 59 nW. Power broad- 
ening corrections were applied to the raw linewidths, but 
these never exceeded 16%. 

Linewidth results for two NV centers in the CVD 
sample are summarized in Fig. |4j For T < 10 K the 
linewidth is fairly constant and not much larger than 
the spontaneous emission lifetime limit of 13 MHz. How- 
ever, for r > lOK the linewidth increases rapidly, fol- 
lowing approximately a dependence. The solid curve 
shows a fit using 7(T) = 70 + C2rT^ up to 100 K with 
70 = 27r X 16.2 MHz, C2 = (9.2 ± 0.5) x 10-^K-^ 
and r = (12.5 ns)~^ as above [24]. Fig. |4] also in- 
cludes spectrometer-based PL linewidth measurements 
obtained from single NV centers, extending the experi- 
mental results to 300 K. Comparing the polarization re- 
laxation and line broadening results, we find that the 



coefficients ci and C2 differ only by a factor of two. In 
these measurements this difference is not statistically sig- 
nificant, however we note that one possible cause for a 
difference could be the different strain environments in 
the two samples. 

To explain this observed dependence, we propose 
a two-phonon Raman scattering process mediated by a 
linear electron-phonon interaction. The Hamiltonian is 
given hy H = He + Hy + H' in which Hg. is the NV elec- 
tronic Hamiltonian and is the vibronic Hamiltonian. 
H' is the electron-phonon interaction and in the weak 
coupling limit can be written in the symmetry adapted 
form (excluding A2 vibrations for brevity) |25j 

k 

+ Y.^k{l^x{ak. + al^) +Uy{aky + aly)) (1) 

k 

in which the first term relates to the interaction with Ai 
vibrations and the second with E vibrations. V^'" gives 
the strength of the coupling to distortions of F symmetry 
for phonon momentum k. Operator a (a^) is the phonon 
annihilation (creation) operator. lA^ are electronic oper- 
ators transforming as row 7 of F symmetry. For strain 
applied along the NV's reflection symmetry axis, in 
the {Ex, Ey} basis have the values [25] 

= ( J ? ) ' = ( J -1 ) ' = " ( ? J ) 

Considering the wavefunctions to first order in the 
electron-vibration interaction, for the case of Ai vibra- 
tions there is no mixing between the E,j. and Ey elec- 
tronic states and as a consequence the polarization se- 
lection rules for the Ai phonon sidebands are the same 
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as for the ZPL. For vibrations of E symmetry, it can 
be seen from Uy that an Ey admixture does oc- 

cur and the emitted polarization can be reversed. Thus, 
the finite polarization visibility observed in Fig. [2}: at 
very low temperatures indicates that there is a significant 
electron-vibration interaction associated with vibrations 
of E symmetry. This is the first evidence of a D JT effect. 

At higher temperatures the terms in the electron- 
vibration interaction associated with E vibrations can 
also induce transitions between the E^, and states 
through a two-phonon Raman process. By using a sym- 
metry adapted approach the rate for a transition Ey 
Ex with an increase in a quantum of y-vibration and de- 
crease of a x-vibration is given by 



Wa:y = 27r / dujn^{n^ + l)pE^{Lo)pEy{i^) 
Jo 

^xxx^xyy ^xyy^yxy 



(2) 



in the limit lE^ ~ Ey \ <^ kT. In Eq.[2] Z^,,^ = {E,\Ui\Ej) 
and the phonon occupation number is given by — 
{exp[hLu/kT] — 1)^^ for phonon energy uj. The phonon 
density of states of the x- and y- vibrations, given by 
Pe-,{uj) — pEy{^), are both proportional to lo^ in the 
Debye approximation. The linear electron-phonon cou- 
pling varies as ^/o; in the absence of resonant phonon 
modes [13]. Thus, the integrand in Wxy varies as 
which leads to a dependence jMj . This is very differ- 
ent from the Raman rate associated with the coupling 
of vibrations to non-degenerate states which varies as 
[27]. The experimentally observed dependence in 
the NV polarization experiment is the second evidence for 
coupling to a degenerate vibrational mode, or a DJT ef- 
fect, in the low strain NV system. When the strain split- 
ting is larger, as with NV3-5, the strain further degrades 
the selection rules and in addition one-phonon relaxation 
results in an asymmetry in the visibility curves for the 
two pump polarizations. This is qualitatively what is 
observed in Fig. [2|:. 

The Raman relaxation shortens the electronic lifetime 
in cither ^E state level resulting in a dephasing of the op- 
tical transition and an increased optical linewidth. The 
observed broadening of the ZPL and the similar mag- 
nitudes of ci and C2 strongly suggest that the DJT effect 
also accounts for the optical dephasing at low T. If the 
linewidth is solely determined by the '^E population re- 
laxation we would expect ci—C2 = 2Wxy/rT^. 

The established presence of a DJT effect in the neg- 
atively charged NV system is significant as the effect is 
known to substantially diminish interactions such as spin- 
orbit, orbital Zeeman, and the response to external per- 
turbations such as stress and electric field [251 12H]- It 
can, therefore, be understood why the effects of many of 
these interactions are found to be small in the NV system. 
The DJT effect may also play a role in the fast popula- 



tion transfer between the excited states observed at room 
temperature. This fast transfer effectively averages the 
two orbital states resulting in a significant simplification 
of the excited state spin resonance structure observed 
at room temperature [55J [301 131] • Finally, by improv- 
ing our basic understanding of the primary decoherence 
mechanism affecting the optical transitions, these results 
should allow much more accurate performance estimates 
for quantum information processing applications includ- 
ing schemes designed to mitigate the effects of excited- 
state dephasing j32j . 
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